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Abstract 

New experimental results, and a plausible theoretical understanding thereof, are presented for 
the flow-induced currents and voltages observed in single-walled carbon nanotube samples. In our 
experiments, the electrical response was found to be sublinear - nearly logarithmic - in the flow 
speed over a wide range, and its direction could be controlled by an electrochemical biasing of the 
nanotubes. These experimental findings are inconsistent with the conventional idea of a streaming 
potential as the efficient cause. Here we present Langevin-equation based treatment of the nanotube 
charge carriers, assumed to be moving in the fluctuating field of ions in the flowing liquid. The 
resulting "Doppler-shifted" force-force correlation, as seen by the charge carriers drifting in the 
nanotube, is shown to give a sublinear response, broadly in agreement with experiments. 
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Single-walled carbon nanotubes (SWNT) in contact with a flowing liquid provide a unique 
microfluidic system that offers a large interfacial area of intimate atomic contact between 
the liquid and the solid substrate. This can lead to a strong coupling of the charge carriers 
in the nanotube to the particles in the flowing liquid, more so if the liquid is polar or ionic 
in character. The effect of this coupling is expected to be further enhanced due to charge 
carrier entrainment because of the quasi-one dimensionality of the conducting nanotubes. 
Recently, the flow of a variety of liquids over SWNT bundles was studied experimentally, 
and was found to generate voltage in the sample along the direction of the flow Quite 
unexpectedly, however, the dependence of the voltage on the flow speed was found to be 
sublinear, and could be fitted to a logarithmic form over five decades of variation of the 
speed. There has been an attempt at explaining away this flow-induced voltage in electroki- 
netic terms as the streaming potential that develops along the flow of an electrolyte through 
a microporous insulator |2]. Earlier Krdl and Shapiro 3] proposed that the shear stress of 
liquid flow transfers momentum to the acoustic phonons of the nanotube and the resulting 
"phonon wind" produces a current of carriers in the nanotube. They also suggested, quali- 
tatively, that the fluctuating Coulomb fields of the ions in the liquid could drag directly the 
carriers in the nanotubes. To produce a flow-rate of 1 cm/s in an inter-nanotube channel of 
20nm width, as demanded by the first mechanism j^] requires an enormous pressure head, 
about 10 13 dynes/cm 2 . The second mechanism jsj gives currents of order femtoAmperes. 
In both mechanisms [sj, the effect is linear in the flow rate, in contrast to the experimental 
findings [1] . The mechanism we offer is related to the second idea of [3] but requires neither 
localization of carriers nor drag at the same speed as the ions. The main aim of this paper 
is to report (i) our new experimental results of the measurements of the the short-circuit 
electrical current as a function of the flow-speed, (ii) ways to control the magnitude and di- 
rectionality of the flow induced current; and (iii) a general theory for the electrical response 
consistent with these new observations. Our measurements clearly rule out the electrokinetic 
mechanism based on the idea of a streaming potential. Our theory is based on a stochastic 
treatment of the nanotube charge-carriers assumed to be moving under the influence of the 
correlated ionic fluctuations which are advected by the liquid flow. It is broadly consistent 
with our experiments. 

Figure ^ shows a schematic sketch of the construction of the flow sensor. SWNT bundles 
prepared by arc discharge method J] are packed between two metal electrodes. The nan- 
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otubes are kept in their place by a supporting insulating substrate. The electrical signal is 
measured along the flow direction (ul) as shown Fig. ^ The other experimental details are 
as in A sensor with a minimal contact resistance of ~25fi (found from four- probe mea- 
surements) was used in the experiments so that the short-circuit current could be measured. 
The short-circuit current (open-circuit voltage) was measured by connecting the microam- 
meter (millivoltmeter) across the SWNT sample. The resistance (two-probe) of the device, 
measured with the sensor dipped in the liquid was found to be ~ 70 fl Figure El shows the 
dependence of the induced voltage and current on the flow velocity u^. The solid line is a fit 
to the empirical relation I = ai log(/3j ul + 1), with ai = 0.02/iA and /3j = 4.8 x 10 4 s/cm. 
The voltage also fits the empirical relation V = ay log(/3y ul + 1), where ay = IA/iV and 
j3y = 4.8 x 10 4 s/cm. It can readily be seen that ay = aj x R, i.e., the resistance encountered 
is precisely the resistance (2 probe) of the device. This is an important point to note: if an 
electrokinetic mechanism were operating, the resistance obtained would have been orders of 
magnitude higher, i.e., equal to that of the electrolyte (~ 0.1MQ) p}. This in itself rules out 
quite decisively the electrokinetic mechanism of voltage generation. Next, we consider the 
measured direction of the flow induced current with respect to the flow direction as a func- 
tion of the bias voltage Vb (see inset of FigJ3J). This potential biases the SWNT with respect 
to the Au-reference electrode immersed in the flow chamber close to sample as shown in the 
inset of Fig.Q. The dependence of the sign and the magnitude of the flow- induced voltage 
on Vb for an aqueous solution of 0.01 M KC1 (conductivity 1.4 mS/m) and for a fixed flow 
speed of 0.04 cm/s is shown in Fig.flHJ). It is seen that the flow-induced signal is positive, 
i.e., I is anti-parallel to ul when Vb is positive, and the sign of the signal is reversed, i.e., 
I is parallel to ul, for Vb negative, consistent with our theory. Further, in our theoretical 
mechanism, based on fluctuations, the carriers in the nanotube drift necessarily in the same 
direction as the flow velocity ul, quite independently of the sign of the ionic charge. Thus, 
for the current / to be parallel (antiparallel) to ul the charge carriers in the nanotubes need 
to be holes (electrons). When the nanotubes are biased positively, the anions (CI - , OH~) 
move closer to the SWNT, localizing holes on the carbon and making electrons available for 

n 

flow-induced current. Similarly, holes are liberated when the bias is negative jjj. As the 
bias voltage is increased the number of carriers participating in the flow-induced current 
will increase as shown in Fig.Q. Thus, the results obtained by electrochemically biasing 
the nanotubes are naturally consistent with our mechanism of (coulombic) forcing of the 
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nanotube charge carriers by the liquid flow. The dependence of the flow-induced signal on 
the concentration of different types of ions in the liquid is, however, found to be complicated 
and non-monotonic. These details will not be addressed here. 

We now turn to our theoretical mechanism for the flow- induced current/voltage. This 
can be understood qualitatively in terms of three physically distinct but related ideas: (a) 
Induced friction: The fluctuating charge density of the ions close to the nanotube couples 
couloumbically to the charge carriers in the nanotube and, therefore, offers a friction to the 
motion of these charge carriers (in addition to the Ohmic friction intrinsic to the carbon 
nanotubes). This, of course, follows directly from the fluctuation-dissipation theorem; (b) 
Flow-induced drag: In virtue of the above frictional coupling, an imposed liquid flow drags 
the charge carriers along through the nanotube; (c) Reduction of induced friction at high flow 
speeds: The space-time correlated couloumbic fluctuations, inherent to the liquid electrolyte, 
are advected by the liquid flow, and thus get Galilean boosted (Doppler shifted) as seen in 
the mean rest frame of the drifting carriers in the nanotube. Correspondingly, as we will see, 
the friction they offer to the motion of the charge carriers in the nanotube diminishes with 
increasing flow speed. This is crucial to the observed sublinear dependence of the charge 
drift- velocity (electrical response) on the liquid flow speed. With the above in mind, we will 
now derive these frictional effects, first from a heuristic argument, and then analytically 
from a Langevin-equation treatment. 

Consider a nanotube placed along the the z-axis with the liquid flowing parallel to it. 
We model the nanotube as a classical one- dimensional (ID) conductor with diffusive charge 
transport. In the steady state, let ul and ud be, respectively, the velocity of liquid(L) flow 
and the induced drift velocity of charge carriers in the nanotubes, all measured relative to 
the nanotube lattice. Now, in the absence of any frictional coupling to the liquid flow, the 
drift velocity u D will tend to relax to zero, i.e., to rest with respect to the lattice, with 
a relaxation time td characteristic of the nanotube resistivity. Similarly, if we were to 
"switch off" the resistive coupling to the lattice, the drift velocity u D would relax to the 
liquid velocity ul, i.e., to the rest frame co-moving with the liquid, because of the frictional 
coupling (drag), with a relaxation time tl- Now, therefore, in the presence of both these 
frictional influences, the drift velocity ud will assume a steady-state value up < ul that 
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satisfies u d /t d = (u L - u d )/t l giving 

u D = u L /(l + r L /r D ). (1) 

Equation JI} is merely a restatement of the condition of frictional force-balance in the steady 
state. It would appear to give an induced short-circuit current (equivalently, an open circuit 
voltage via the nanotube resistance) along the nanotube which is linear in the flow velocity. 
The nonlinearity is, however, really hidden in the ul dependence of the relaxation time 
Tl that we will now try to make explicit. It may be noted here that we are assuming, 
for simplicity, a uniform liquid flow without the hydrodynamic complications of a no-slip 
boundary condition. 

In a simple caricature of the real situation then, consider the ionic density in the liquid, 
fluctuating thermally and flowing past the nanotube at a mean velocity u^z,, producing 
thereby a fluctuating couloumbic potential </>(r, t), at a point r at time t. We are, of course, 
interested in the case of r lying on the z axis i.e., r = (0,0, z) (in the ID nanotube). 
For the space-time correlation function (0(0, O)0(r, £)) = Go(r,t) in the mean rest frame 
of the ions, the charge carriers in the nanotube see this correlation Galilean boosted to 
G(r,t) = Go(r — zvt,t) with v = ul — ud- This Galilean boost (Doppler shift) is the key 
physical point of our treatment. At Ul — 0, the fluctuation-dissipation theorem (FDT) tells 
us that the coefficient of the zero-frequency friction to the motion of the charge carriers 
in the nanotubes, arising from the ionic thermal fluctuations, is proportional to the time 
integral of this on-site force-force correlation function. If we assume this relation even for 
«i ^ we have 

/oo 
(eE z {r-zvt,t)eE z {r,0))dt. (2) 
-oo 

Here E z is the z-component of the coulombic (electric ) field due to the ions; m e is the 
mass of the charge carrier with e the electronic charge; ks is the Boltzman constant, and T 
the absolute temperature. We re- write the right hand side of Eqn.(J2J) in Fourier (q)- space, 
expressing the above force-force correlator in terms of the ionic charge-densities p(r, t) using 
E z (q,t) = —iq z (f>(c[,t) and — g 2 0(q) = ep(q, t)/e, where e is the solvent dielectric constant, 
and obtain straightforwardly 
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where x = v/Dk, and po is the mean ionic number density. In Eqns. J3J) and (J3J we 
have used the Debye -screened form for the static charge structure factor Sq = (p q p- q ) = 
q 2 j (q 2 + k 2 ) with screening length k _1 as the inverse of the Debye screening length and a 
diffusive form l/r q = Dq 2 |7| with D the ionic diffusion constant. It can be seen at once 
that Eqn.(@J) , taken in conjunction with Eqn.(j2J), gives a drift velocity U£,(and therefore the 
short-circuit current) as a sublinear function of the flow velocity ul- This sublinearity is a 
generic feature of this mechanism, and is clearly seen in inset of Fig.© in which we have 
plotted the induced current versus the liquid flow speed derived directly from Eqn.(HJ) for 
certain choice of parameters. Rather than attempting to justify this use of the regression 
formula in a nonequilibrium context, we turn instead to a direct evaluation of the drift speed 
starting from a Langevin equation. We recover the same sublinear functional form for the 
response as in Eqn.(JSJ). 

The one-dimensional (ID) position z(t) of a tagged charge carrier in the nanotube obeys 
the overdamped Langevin equation 



where F(z(t),t) is the fluctuating force the ions exert on the carrier at z(t) in the nanotube, 
and / and ( are the thermal Gaussian white noise and friction intrinsic to the nanotube, 
with (f(0)f(t)) = 2kBT(5(t). Strictly speaking, F(z(t),t) should be a function of the 
instantaneous separations between the positions z(t) of the tagged carriers in the nanotube 
and those of all the ions in the ambient liquid. We would like to replace this complex many- 
body problem by an effective Langevin equation for the dissipative motion of the charge 
carriers in the nanotube, without invoking the FDT which, strictly speaking, holds only in 
the absence of the flow. Our treatment nonetheless captures the qualitative physics of the 
dissipative entrainment of carriers by the moving ions. For this, we simply take F(z, t) to 
be a Gaussian random force field, with zero mean and a two-point correlation C(z, t). Then 
from Eqn.(jHJ), the drift velocity of the carrier Ud = (dz(t)/dt) = (l/()(F(z(t),t)) can be 
rewritten, using Novikov's theorem Q in relation to the Gaussian noise F(*,t), as 



Using Eqn.(jSJ) to evaluate the functional derivative in Eqn.(jHl), and writing C(z,t) in terms 



C =F(z(t),t) + f(t) 



(5) 
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of its Fourier transform C q (t) yields 



u D 



dq 



dt{e iq ^- z ^)iq z C q {t) (7) 



C 2 J-a (2tt)3 Jo 

where A is an ultraviolet cut-off of the order of an inverse ionic diameter. As in the 
preceding heuristic treatment, let us take the correlation C q (t) of the ions to be the Galilean 
boost, with velocity ul, of an equilibrium correlation function C q (t), with a relaxation time 



C q (t) = C q {t)e~ iq * ULt = C° q e- iqzULt e- t/r ", (8) 

where C q is the equilibrium equal-time correlation function of the force fluctuations. This 
form, despite its undeniable limitations, is the simplest way to capture the basic physics 
of ions moving past the nanotube, and admits an essentially analytical treatment. As 
before, the force-force correlation C q can be expressed in terms of the ionic charge-density 
correlation, which is known as an input from the liquid state (dilute ionic solution) theory, 
namely that C q oc (q 2 z /q A )S q with the ionic charge structure factor S q = q 2 /(q 2 + k 2 ). 
Note the factor (ql/q^) arising from the gradient (d/dz) and the Laplacian(V 2 ) in q- space. 
Replacing z(t) in Eqn.(J7J) by its mean ur>t for simplicity, we obtain the compact expression 

«* = «« t. 7^(3^) {^rX 0» 



-a (27r) 3 \q 2 z v 2 + T q 2 J \q 2 + k 

for the drift velocity of the charge carriers in the nanotube, where as before v = ul — ud, a 
is a lumped constant of proportionality that depends on the parameters of the liquid-state 
correlation function input used above. With the ultraviolet cut-off (A) set to infinity, and []| 
with \jr q — Dq 2 -, Eqn. Q has precisely the form of Eqn.(|T|) taken in conjunction with the 
Eqn.(|Hl. allowing us thereby to identify the integral on its right-hand side essentially with 
t d /t l . This gives us an expression for the flow-speed dependence of t l , and thus finally an 
analytic expression for the charge drift velocity («d) as a function of the liquid flow velocity 
{ul). This reaffirms our heuristic argument given at the beginning. 

We close by summarizing the main points of our work. First, on the experimental side, 
we have clearly shown that the liquid flow produces not only a voltage (i.e., not merely a 
capacitive charging), but a short-circuit current as well in the nanotube; that both have a 
sublinear dependence on the imposed flow-speed; and that the voltage/current ratio cor- 
responds to the nanotube sample resistance. These observations are incompatible with an 
electrokinetic origin for the (electronic) current in the nanotube. On the theoretical side, 



we have proposed a theory wherein the current is essentially a statistical consequence of the 
flow-induced asymmetry in the correlation of the ions, in the ambient fluid as seen by the 
charge carriers in the nanotube. Importantly, our theory predicts in general a sublinear be- 
havior for the electrical response, with a linear regime at only the smallest values of imposed 
flow. The extended logarithmic regime seen in experiments can presumably be rationalized 
in detail with particular forms for the correlation function (C°) and the relaxation time 
(r g ), as inputs to be taken from the liquid state (dilute ionic solution) theory. Moreover, 
a realistic treatment will require taking into account details of the complex, hydrophobic, 
inter-nanotube micro-fluidic environment of our mat samples. Thus, very specifically, the 
no-slip boundary condition would imply a decreasing velocity of the flow nearer the nanotube 
(the shear flow). This decrease in the flow velocity will,however,be offset by the corresponding 
increase in its effectiveness (via the screened coulombic forcing) closer to the interface. The 
resulting levelling is expected to broaden the sublinear response and thus improve agreement 
with the experiment. Our main point, however, is that the experiments show sublinear be- 
havior, which is inconsistent within existing theories of flow-induced voltages and currents 
in nanotubes, and that our approach naturally and inevitably leads to strong sublinearity. 
Finally, we emphasize that the flow-induced asymmetry of the random fluctuations is key to 
the charge- carrier drift (drag) mechanism in our theory. In this broad sense our approach 
here subsumes the asymmetric fluctuating ratchets invoked earlier in a general way. In 
this connection, reference must be made to the idea of a drag (shear) induced by the relative 
motion between material surfaces, where the Doppler shifted and aberrated photonic fluc- 
tuations, e.g. zero-point photons, have been invoked very effectively 0, • We would also 
like here to add that more than one mechanisms could very well be at work in these sys- 
tems. For example, in a recent publication Persson et al 11J have invoked a combination of 
frictional stick-slip and barrier-hopping to explain the observed phenomenon of flow-induced 
voltages in SWNT. The one we propose here seems particularly robust and general, and we 
look forward to experimental tests, especially of the predicted saturation of the electrical 
response at high flow speeds. 
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FIG. 1: Schematic sketch of the nanotube flow sensor placed along the flow direction (ul). SWNT 
bundles are packed between two metal electrodes. The insulating substrate keeps the SWNT in 
place. The electrical leads are taken out from the metal electrodes. 



10 




(cm/s) 

FIG. 2: Voltage (full circles) and current (open circles) as functions of flow speed ul- The solid 
line is a fit to the logarithmic function as explained in the text. Inset shows the theoretical 
plot of current (/ = neu^A) versus flow speed based on Ean. (Jll4j) for typical choice of parameters: 
T=300K; e=80( CGS units); Dk = 10~ 4 cm/s; td = 10 -16 s; po = 10 13 cm -3 ; charge carrier density 
in nanotubes (n = 10 18 cm" 3 ); cross-sectional-area (A) = 10 -3 cm 2 . The strong sublinearity is 
clearly seen. 
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FIG. 3: Flow- induced voltage as a function of bias Vb- Inset: schematic of electrochemical biasing 
of the nanotubes; CE is the counterelectrode.The solid line is a guide to the eye 
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